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ABSTRACT 



We constrain the velocity power spectrum shape parameter T in linear 
theory using the nine bulk-flow and shear moments estimated from four re- 
cent peculiar velocity surveys. For each survey, a likelihood function for T 
was found after marginalizing over the power spectrum amplitude crg,^l m 6 
using constraints obtained from comparisons between redshift surveys and 
peculiar velocity data. In order to maximize the accuracy of our analyses, the 

■ velocity noise ex* was estimated directly for each survey. A statistical analysis 

of the differences between the values of the moments estimated from different 
surveys showed consistency with theoretical predictions, suggesting that all 

\ the surveys investigated reflect the same large scale flows. The peculiar ve- 

■ locity surveys were combined into a composite survey yielding the constraint 
£ — | T = 0.13i ; 05 . This value is lower than, but consistent with, values obtained 

■ using redshift surveys and CMB data. 



Subject headings: cosmology: distance scales - cosmology: large scale structure of the universe - cosmology: obser- 
vation - cosmology: theory - galaxies: kinematics and dynamics - galaxies: statistics 
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1 INTRODUCTION 

Observations of the large scale peculiar velocity field provide an important tool for probing mass fluctuations on 
~ 100 h _1 Mpc scales (h is the Hubble constant in units of 100 km s" 1 Mpc -1 ). Analyses of peculiar velocity 
surveys used to constrain the amplitude of mass power spectrum (Freudling et al. 1999; Zaroubi et al. 2001) are 
complementary to those that employ redshift surveys alone (e.g. (Percival et al. 2001; Tegmark Hamilton & Xu 
2002)) or in combination with CMB data (Elgaroy, Gramann & Lahav 2002; Sanchez etal 2006). Although initially 
peculiar velocity surveys using redshift independent distance indicators were sparse and shallow (Faber & Jackson 
1976; Tully & Fisher 1977) , large, homogeneous redshift-distance samples of galaxies and clusters have become 
increasingly common. Early analyses of redshift-distance surveys (Aaronson et al. 1982; Lynden-Bell et al. 1988) 
led to the development of powerful statistical methods for the analysis of peculiar velocity data (Dressier etal. 1987; 
Kaiser 1988; Feldman & Watkins 1994; Strauss & Willick 1995; Watkins & Feldman 1995), but were hindered by 
shallow and sparse samples. However, with the advent of larger and better samples (Feldman et al, 2003; Pike & 
Hudson 2005; Park & Park 2006; Sarkar, Feldman & Watkins 2007) it has become increasingly clear that peculiar 
velocity catalogues can play an important role in the determination of cosmological parameters. 

Recently, several studies have compared peculiar velocity surveys directly to the density field as determined 
from redshift surveys and put strong constraints on the combination of parameters ogA^ 6 (for a summary, see 
Pike & Hudson 2005). These studies have involved peculiar velocity data derived from several different distance 
measures, as well as three different redshift surveys, and yet have produced remarkably consistent results. When 
one considers that the constraints from these studies are also consistent with those derived from other types of 
studies, it seems reasonable to consider our knowledge of agQ^ to be firmly established. 



© 0000 RAS 



2 Watkins & Feldman 



The combination of parameters o&Sl^ essentially determines the amplitude of the velocity power spectrum. 
The shape of the power spectrum can also be constrained using peculiar velocity data, but efforts of this type have 
typically tried to constrain the amplitude and shape together, resulting in constraints on a two-dimensional space 
of parameters (see, e.g. Jaffe & Kaiser 1995, Borgani et al. 2000). In this paper we use peculiar velocity surveys to 
put constraints on power spectrum parameters directly through linear theory. By using the constraint on crs^m 6 as 
a prior, we are able to calculate the likelihood function for the "shape" parameter, F, of the power spectrum by 
itself, resulting in a simple one-dimensional constraint on this parameter. 

Given that our analysis uses linear theory, it is important that we properly filter out small-scale flows due 
to nonlinear effects that are present in peculiar velocity surveys. For this reason, we follow Jaffe & Kaiser (1995) 
in analyzing surveys using only the two lowest order moments of the velocity field, i.e. the nine bulk flow and 
shear moments. For the sizes and depths of the surveys that we are considering, the bulk flow and shear moments 
probe velocity modes with wave numbers well above the scales where nonlinear effects are thought to be important. 
Significantly, the bulk flow and shear moments are also the velocity moments that have the highest signal-to-noise 
and thus can be determined most accurately; indeed, little useful information is lost by discarding higher order 
moments. 

We apply our analysis to an extensive set of peculiar velocity surveys. Each of these surveys employs a different 
method of distance estimation, and target different populations of galaxies. We expect that each survey is affected 
by nonlinearities in different ways; this could be reflected as differing amounts of small-scale motion superimposed 
upon the large-scale linear flow reflected in the bulk flow and shear. These small-scale motions can be quantified by 
a standard deviation a* of the velocities remaining after the bulk flow and shear have been subtracted out. While 
a, is typically given a fixed value of around 300 km/sec in this type of analysis (Kaiser 1991; Feldman & Watkins 
1994), in order to improve the accuracy of our study we have calculated the maximum likelihood value of a* directly 
for each survey using an iterative method. 

An important question, then, is how well does the information about large-scale flows contained in these 
different surveys agree? In order to answer this question quantitatively, we calculate the covariance matrix for the 
differences between the values of the bulk flow and shear moments for two surveys. The covariance matrix, together 
with the measured differences, allows us to calculate a % 2 that reflects the probability that both surveys reflect the 
same underlying large-scale flow. This probability is most useful when surveys are similar in their characteristics, 
as the surveys that we consider are; two surveys that probe the velocity field in different ways can have a high 
probability of agreement even when the values of their bulk flow and shear moments are quite different. 

In § 2, we give the details of our likelihood analysis. In § 3 we discuss the power spectrum parametrization we 
use. In § 4, we describe the peculiar velocity surveys used in the analysis. In § 5 we present our results, and in § 6 
we discuss them and compare them to constraints on the shape parameter derived from other types of data. 



2 LIKELIHOOD ANALYSIS 

Small scale motions of galaxies reflect nonlinear evolution and can depend strongly on the specific population of 
galaxies that are sampled in a survey. However, by forming velocity moments out of weighted sums of individual 
velocities that reflect only large scale motions, we expect small scale motions to average out so that linear theory 
is applicable. In this study, we focus on the measures of the large scale flow given by the first and second order 
moments of a taylor expansion of the velocity field, 

Vi(r) = Ui +Pijr$ H , (1) 

where u is the bulk flow vector and pij is the shear tensor, which can be taken to be symmetric. While the 
interpretation of these measures depend on the specific distribution of galaxies in a survey as well as measurement 
errors, the surveys we consider are similar enough that they form a good basis for comparison (Sarkar, Feldman & 
Watkins 2007). 

The samples that we consider consist of a set of N galaxies, each with a position vector r n and a measured 
line-of-sight peculiar velocity S n = v-f n with individual measurement error a n . Following (Kaiser 1991) , we group 
the bulk flow and shear components into a 9 component vector a p , so that the galaxy velocities can be modeled as 

S n = a p g p (r n ) + 5 n , (2) 

where the 9 component vector g p (r) = (r x ,f y ,f z ,rr x r x , rr x ry,rr x f z ,rr y ry,rfyf z ,rr z r z ). In what follows we take 
our coordinate axes to correspond to galactic coordinates. The deviation from the model, <5„, contains contributions 
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from both small scale motions and measurement errors. We assume that the 8 n values are gaussian distributed and 
have a variance given by a'i + o\ (Kaiser 1988). Here <t» represents the small-scale linear and nonlinear motions 
that are not accounted for by the bulk flow and shear. We shall refer to cr» as the velocity noise. 
Under these assumptions, the likelihood function for the moments is 

For a given survey, we would like to find the maximum likelihood values for the moments a p and for cr*. In 
previous analyses of this type, the value of cr* has sometimes been determined by an educated guess (see, e.g. 
(Kaiser 1988; Feldman & Watkins 1994; Jaffe & Kaiser 1995; Pike & Hudson 2005)). However, we note that since 
surveys generally sample a specific population of galaxies, each of which can be differently affected by nonlinear 
overdensities, we expect that each survey has, in principle, a different value for cr*. Thus, here we determine the 
value of a* directly from each survey. We find the maximum likelihood values for the moments a p and cr* using an 
iterative method. We start by making a "guess" of cr, = 300 km/sec. Treating cr* for the moment as constant, the 
maximum likelihood solution for the a p is given by 

A -i \ - g q (rn)S n 

c n "T u * 

n 

where 

A _ \ " ffp(rn)ffq(r n ) 

n 

Using these estimates in Eq. (3) and treating them as constant, we then find the maximum likelihood value of 
cr* , which now replaces our original guess and can be used to calculate a refined set of maximum likelihood values 
for the a p . This process is repeated until the estimates converge, which in practice only takes two or three iterations. 

In order to compare our estimates of the moments a p to the expectations of theoretical models, we calculate 
the covariance matrix, which from Eq. (4) can be written as 

R pq = (a p a q ) — A pl A qs „ 2W_2 T _2N (SnSm), (6) 

\"n ~T" °* ) m ~r~ (J * ) 

n.vri 

where {S n S m } can be written in terms of the linear velocity field v(r) and the variance of the scatter about it, 
{S„S m ) = (?„ • v(r„) f m ■ v(r m )) + 8 nm (al + al). (7) 
Plugging Eq. (7) into Eq. (6), the covariance matrix reduces to two terms, 

R p q ~ R p q ~\~ R p q • (8) 

The second term, called the "noise" term, can be shown to be 

R p q ^-pq ' (^) 

The first term is given as an integral over the matter fluctuation power spectrum, P(k), 

q1.2 roo 

R ( v1 = TTT dk W 2 P q(k)P(k), (10) 



2tt 2 Ju 

where the angle-averaged tensor window function is 



Given a peculiar velocity survey and the values of its 9 moments, we can write the likelihood of a theoretical 
model used to calculate the covariance matrix as 

C = j^jTTJ ex P (-\ a p R pq la i) ■ ( 12 ) 

As we describe below, we use this equation in order to place constraints on the parameters of cosmological models; 
in particular T, the parameter that determines the shape of the power spectrum. 

Both the bulk flow and shear moments probe primarily large-scale motions and thus have window functions 
that are peaked at large scales. However, our prior constraint on agQ 1 ^ 6 essentially fixes the amplitude of the power 
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spectrum on relatively small scales. Since the shape parameter Y controls the relative distribution of power between 
large and small scales (smaller V corresponds to relatively more power on large scales), increasing F generally results 
in smaller values for the elements R pq of the covariance matrix. Thus Eq. (12) shows that surveys with large values 
of the bulk flow and shear moments generally favor smaller values of T. 

The window function for a given survey carries information about how the moments of that survey sample 
the power spectrum. Since each survey has a unique window function, the values of the moments are not strictly 
comparable between surveys (Watkins & Feldman 1995). However, since the volumes occupied by the surveys we 
consider overlap strongly, we expect the values of the moments of these surveys to be highly correlated. In order to 
quantify the agreement between two different surveys, say, survey A and survey B, we use the covariance matrix 
for the difference between the values a p and a p of the moments for the two surveys, 

r>A — B // A B\/ A B\\ ryA , nB U AB H AB 

tXpq \\^p ^p )\dq )/ ^pq < *^pq **>pq ^qp , V / 

where the cross-terms are given by 



R tf = -^J o dk (W AB ) 2 pq (k)P(k), (14) 



and 



,ABs2 lu s _ , ,A r l, ,fl r lV T 3i( r n)ffs( r m) 



(W )pq(k) - (A ) pl (A )„ ^ + {a ^ 2){{a £ )2 + fl 

x j ( f ™ ' ^ r„ • k) exp [ik ■ (r^ - r„)] . 

(15) 

Note that we have assumed here that the nonlinear contributions to the velocities of the galaxies represented 
by cr» in the two surveys are uncorrelated. This is not likely to be true in reality, since galaxies in the same 
local neighborhood are affected by the same small-scale flows. However, this will always cause us to underestimate 
the expected amount of correlation between surveys. Thus our results on how well two surveys agree should be 
considered as upper bounds. 



3 POWER SPECTRUM PARAMETERS 

As shown in Eqs. (10) and (14), in linear theory the variances for the velocity components a p are given by integrals 
over the power spectrum P(k) multiplied by the factor fi^; 2 . The power spectrum P(k) can itself be modeled as 
an initial power law k n , where we assume the usual n w 1, times the square of the transfer function T(k), so that 
P(k) oc kT 2 (k). We set the overall amplitude of the power spectrum in the usual way through the constant ag, the 
amplitude of matter density perturbations on the scale of 8 /i _1 Mpc. The transfer function is generally parametrized 
in terms of the "shape" parameter T ^ fl m h. We follow Eisenstein and Hu (1998) in writing the transfer function 
as 

T{k) = L + C (fc/F)2' 

L (k) = ln(2e+ 1.8(fc/r), (16) 
7^11 

CJk) = 14.2 + ; - , 

y ' 1 + 62.5(fc/r) ' 

which is valid in the zero-baryon approximation. While Eisenstein and Hu (1998) include a more accurate model 
for the power spectrum that includes the effects of baryons, this model depends on several more parameters and 
does not have a constant V. We have found that there is not enough information in the peculiar velocity data to 
break the degeneracies between these parameters, and that a two parameter model is sufficient given the precision 
of our results. We will revisit the issue of baryons in the discussion section (§ 6) when we interpret our results. 



4 PECULIAR VELOCITY SURVEYS 

We have applied our analysis to four peculiar velocity catalogues. These catalogues vary in their sample size, depth, 
distance measurement method, and typical measurement errors. Our subset of the spiral field /-band (SFI) survey 
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(Haynes et al. 1999a, 1999b) consists of 1016 galaxies at distances d < 60 7i _1 Mpc with velocities |v| < 2000 
km/sec. The distance measurements, obtained using the /-band Tully-Fisher (TF) relation, have been corrected for 
Malmquist and other biases as described in Giovanelli et al. (1997). The distance cut was made to avoid the bias 
toward inwardly moving galaxies near the redshift limit of the survey. The cut on large velocity galaxies was made 
to prevent them from having an undue effect on our results. We have tested that our results are not sensitive to 
the precise way in which these cuts are made. Distance errors on this sample are of order 16%. 

Our subset of the Nearby Early-type Galaxy (ENEAR) survey (da Costa et al. 2000a, 2000b) contains 535 
galaxies at distances d < 70 h _1 Mpc with velocities \v\ < 2400 km/sec. The D n — a distance estimates in this 
catalogue have been corrected for inhomogeneous Malmquist Bias (Bernardi et al. 2002) . A correction has also been 
applied for the bias toward inwardly moving galaxies at the redshift limit of the catalogue; however, given the large 
size of this correction for galaxies near the limit, we chose to apply a distance cut in order to avoid the effects of 
this bias and its correction on our sample. Distance errors for this sample are of the order of 18% for individual 
galaxies with groups of size N having their distance errors reduced by a factor of ^/N. 

Our subset of the /-band surface brightness fluctuation survey of Galaxy distances (SBF; (Tonry et al. 2001)) 
contains 257 E, SO, and early-type spiral galaxies out to a redshift of about 4000 km/sec. We have removed galaxies 
with percentage distance errors (> 20%) as well as three high-velocity galaxies with redshifts > 4000 km/sec. More 
than half of the galaxies in our sample have distance errors below 10%, with the smallest errors being less than 5%. 
With errors of this magnitude, our sample is immune from the effects of Malmquist bias (Tonry et al. 2000) . 

Finally, we have included a small sample of 73 supernovae of Type la (SNIa; Tonry et al. 2003) whose redshifts 
extend to 8000 km/sec. Typical Distance errors for this sample are about 7%; one supernova whose distance error 
exceeded 20% was not included. As with SBF, this sample should be immune from Malmquist Bias. 



5 RESULTS 

The first step in our analysis is to obtain estimates for the moments a p and the velocity noise <r* for each catalogue 
by maximizing the likelihood function (Eq. 12) as described above. The calculated moments and cr* for each survey 
are given in Table 1. Note that the estimates for the bulk flow vectors are quite similar, as one would expect for 
surveys that occupy similar volumes. The shear moments are somewhat less consistent due to both the facts that 
these moments are estimated less accurately and also depend more strongly on the details of the sample geometries. 
The estimates of er» are all of the expected magnitude. We shall leave a more detailed examination of the a* 
estimates for the discussion section. 

Table 1: Maximum likelihood estimates for the values of the three bulk flow and six shear moments as described in 
Section 2 (equation 2), and the velocity noise ex*, for each of the surveys of interest and the composite survey. 



Survey Bulk Flow Shear <t» 





1 


km s 


x ) 






(km s 


^pc- 1 ) 






(kms- 1 ) 


SFI 


69.8 


-183 


47.6 


1.35 


-2.94 


3.37 


-3.08 


-3.2 


-7 


413 


ENEAR 


142 


-204 


47 


3.65 


-2.39 


5.11 


-4.23 


-4.3 


-3.86 


386 


SBF 


120 


-309 


206 


9.88 


-3.45 


-2.99 


3.41 


-1.99 


-16 


304 


SNIa 


57.6 


-419 


40.7 


3.39 


-3.25 


5.19 


2.62 


-0.848 


-3.14 


327 


Composite 


98.4 


-216 


90.6 


2.6 


-3.27 


3.58 


-2.95 


-3.84 


-6.44 





Next, we use Eq. (12) to calculate likelihoods for model parameters for each catalogue. From Eq. (16), we 
can see that in the context of our model, the theoretical covariance matrix for the velocity moments is completely 
specified by two parameters; the amplitude, given by os^m , and the shape parameter Y. As was mentioned above, 
the amplitude parameter has been strongly constrained; in particular, by comparisons of peculiar velocity data 
and redshift surveys. Pike and Hudson (2005) summarize these constraints and combine several of them to obtain 
as(f2m/0.3) 0,6 = 0.85 ± 0.05. However, we feel that the averaging done on correlated data is not entirely justified, 
and that the constraint (T8(fi m /0.3) ' 6 = 0.84 ± 0.1, or agf^ 6 = 0.41 ± 0.05, which coincides with that obtained by 
Zaroubi et al. (2002), is more representative of the strength of the constraint that can be placed on the amplitude 
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Figure 1. The likelihood functions for Y obtained from each of the surveys. The surveys give consistent results. The thick 
solid line is the likelihood function for the composite survey. 



using peculiar velocity and redshift data. This constraint is also consistent with measurements of this combination of 
parameters using other methods (for a discussion see Pike and Hudson (2005)). Rather than treating the amplitude 
as a free parameter, we instead chose to adopt this constraint as a prior and to marginalize over it. Thus, given 
a peculiar velocity catalogue, we are able to calculate the likelihood function for the single parameter V which 
determines the shape of the power spectrum. Note that although the amplitude and shape of the power spectrum 
are theoretically related through a common dependence on tt m ; rather than presuppose this dependence we choose 
to treat them as independent parameters. We revisit the interdependence of these parameters in the discussion 
section. 

In Figure 1 we plot the likelihood functions obtained from each of the surveys. These likelihood functions are 
asymmetric and have nongaussian tails. In Table 2 we give the maximum likelihood values of F for each survey 
together with the region around the maximum that contains 68% of the probability under the curve. We also list 
the x 2 a t the maximum likelihood, where 

X = 2J a pRp, q a q- (17) 
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These values show that the peculiar velocity data is quite consistent with the power spectrum model we are 
considering. 

Table 2: The maximum likelihood value of 
r and the \ 2 f° r 9 degrees of freedom for 
each of the surveys. The last entry is the 
results for the composite survey. 



Survey T \ 2 



SFI 


16+ 013 
iJ.ro_o.o7 


10.70 


ENEAR 


u - id -0.07 


10.03 


SBF 


o.ir_°o^ 


7.42 


SNIa 


o.i5_° ;?? 


7.62 


Composite 


O-lStoTs 


10.41 



The fact that all of the surveys that we have considered yield consistent results for the value of F does 
not necessarily indicate consistency in the actual values of their moments. In order to check for more detailed 
compatibility between the surveys we consider the question of whether the differences between the values of the 
individual moments of any two surveys, dp — a B , are consistent with that predicted by the theoretical models, i.e. 
are the measurement errors, the velocity noise, and the differences in how each survey probes the power spectrum 
large enough to explain the differences in the moments. Again, we use a \ 2 analysis; calculating the covariance 
matrix Rp q ~ B as described above we form 

X = ~ a v)\ R vi ) [a q -a q ). (18) 

p,1 

It turns out that the \ 2 calculated in this way does not depend very strongly on F in the region of interest. For 
simplicity, then, we report \ 2 values calculated for the single value of F = 0.14. Other values of F give similar 
results. 

The results of this analysis are tabulated in Table 3. They show good consistency between the catalogues for the 
favored range of F values. Thus the velocity moments of the surveys that we consider agree not only in magnitude, 
but also in value, inasmuch as they are expected to given measurement errors, velocity noise, and differences in 
survey volumes. 



Table 3: \ 2 f° r 9 degrees of free- 
dom for the differences between 
the surveys. 



Surveys 


x 2 


SFI-ENEAR 


4.606 


SFI-SBF 


7.430 


SFI-SNIa 


6.104 


ENEAR-SBF 


8.306 


ENEAR-SNIa 


3.299 


SBF-SNIa 


5.194 



Given our result that the four surveys that we have studied are consistent with each other, it seems reasonable 
to combine them into a composite survey which can then be used to obtain the strongest possible constraint on F. 
Since the different values of <r* for the various surveys reflect differences in the populations and distance measures 
between the surveys, we assign each galaxy in the composite survey the value of <r* of its parent survey. In Figure 
1, we show the likelihood function for F resulting from the composite surveys, with the maximum likelihood value 
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being 0.13±g;g§ . In Table 1 we give the maximum likelihood values for the bulk flow and shear moments for the 
composite survey. In Table 2 we present the maximum likelihood value of V for the composite survey together with 



6 DISCUSSION 

Our results show consistency between catalogues containing galaxies of different morphologies and using different 
methods for determining velocities, confirming that each of these catalogues trace the same large-scale velocity 
field within uncertainties. While previous studies have shown consistency in the bulk flow vectors calculated from 
different surveys (Hudson et al. 1999; Sarkar, Feldman & Watkins 2007), our result is the first to directly compare 
both bulk flow and shear moments. 

The samples we considered do show differences on small scales; in particular, the velocity dispersion about the 
bulk flow and shear motions, represented by cr*, shows a range of values. While these differences in <r* could arise 
from how different galaxy populations respond to nonlinearities, underestimates of measurement errors are also a 
likely source for u». 

Although we have used a two-parameter model of the power spectrum that is strictly valid only for the zero- 
baryon case, where theoretically F = fi m /i, in interpreting our results for F it is possible to include the effects of 
baryons to a first approximation. In particular, Sugiyama (1995) has determined that T scales with baryonic density 
fib as 



The parameters in this formula are tightly constrained by microwave background studies (Spergel et al. 2006); 
specifically, h = 0.732±°°g, Q m = 0.241 ± 0.034, and Q b h 2 = 0.0223lg;^|, which can be combined to give tt b = 
0.041 6±0.0049. Plugging these values into Eq. (19) and propagating uncertainties gives the result F = 0.137±0.025, 
which is clearly consistent with our results. While this model is not as accurate as that of Eisenstein & Hu (1998), 
the latter introduces complications in interpretation due to its T having k dependence. We have found that the use 
of the more complicated model in our analysis did not change our results significantly given the precision that can 
be achieved using the available data. 

The constraint that we have obtained is also consistent with other measurements of the power spectrum. In 
particular, an analysis of the 2dF galaxy redshift survey has found that T — 0.168 ±0.016 (Cole et al. 2005) for blue 
galaxies, whereas in earlier measurement (Percival et al. 2001) they found T = 0.20 ± 0.03. The SDSS collaboration 
found T = 0.213 ± 0.023 for all galaxies (Tegmark et al. 2004) and T = 0.207 ± 0.030 for luminous red galaxies 
(LRG) (Pope et al. 2004). Using CMB data from the WMAP experiment together with the SDSS FRG data set 
Hiitsi (2006) found T = 0.202±°'°|}. 

We have used as a prior the constraint on agflm 6 derived from comparisons of peculiar velocity and redshift 
surveys, which is appropriate for an analysis of velocity data. However, other measurements of this quantity exist, 
in particular from microwave background studies (Spergel et al. 2006). Using the value of fi m from above, together 
with the WMAP measurement of cr 8 = 0.776tJ5;°32 , and propagating errors, we find that cr 8 ^ 6 = 0.330lEJ;°42, 
which is about one standard deviation below the central value we have used in this paper. In order to understand 
the effect of the prior for os^m 6 , we show in Fig. 2 the likelihood contours for the composite survey over the full 2- 
dimensional parameter space of F and agQ^ 6 ■ From the plot it is clear that using a lower the value of us^m 6 would 
lead to a maximum likelihood value for F that is smaller, although not significantly so. The plot also illustrates 
the substantial degeneracy between the two parameters; the ridge that extends diagonally across contains a large 
region where the likelihood is close to its maximum. Without applying a prior, as we have done, it is difficult to 
obtain meaningful constraints from the full 2-dimensional analysis. 

It is noteworthy that combining all of the survey data into a single survey did not improve the constraint 
on T as much as one might have expected. While combining the data did improve the measurements of the bulk 
flow and shear, cosmic variance limits how strong a constraint that even perfect knowledge of these nine quantities 
measured for a single volume can provide. It is important to remember that the power spectrum determines only 
the variances of these nine quantities; thus our situation is similar to that of trying to constrain the variance of a 
statistical distribution from nine numbers drawn from that distribution. In light of this, significant improvements 
of our constraint will be difficult to obtain. Increasing the number of moments by expanding the velocity field to 
higher order is likely to move the analysis into a regime where linear theory is no longer applicable. Deeper surveys, 
whose bulk flow and shear are sensitive to larger scale power, could strengthen our constraint; however, given 



its associated \ 2 - 




(19) 
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Figure 2. The 2-dimensional likelihood contours for the composite survey. The dot indicates the maximum likelihood. The 
square represents the maximum likelihood value for T when marginalizing over agO,^ as described in the text. The triangle 
shows the central values from WMAP as shown in the text. The contour levels contain 68%, 95% and 99% of the total 
likelihood in the region of parameter space shown in the plot. Note that the location of these contours would change if a 
different region of parameter space was selected. 

that measurement errors typically increase linearly with distance, the number of galaxies required for a reasonable 
analysis is prohibitive. 
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